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LINKED EVOLUTION OF GAS AND STAR FORMATION IN GALAXIES OVER COSMIC HISTORY 
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ABSTRACT 

We compare the cosmic evolution of star formation rates in galaxies with that of their neutral hydrogen densi- 
ties. We highlight the need for neutral hydrogen to be continually replenished from a reservoir of ionized gas to 
maintain the observed star formation rates in galaxies. Hydrodynamic simulations indicate that the replenishment 
may occur naturally through gas infall, although measured rates of gas infall in nearby galaxies are insufficient to 
match consumption. We identify an alternative mechanism for this replenishment, associated with expanding su- 
pershells within galaxies. Pre-existing ionized gas can cool and recombine efficiently in the walls of supershells, 
molecular gas can form in situ in shell walls, and shells can compress pre-existing molecular clouds to trigger 
collapse and star formation. We show that this mechanism provides replenishment rates sufficient to maintain 
both the observed HI mass density and the inferred molecular gas mass density over the redshift range < z < 5. 



Subject headings: galaxies: evolution — galaxies: formation - 
general — supernovae: general 



galaxies: ISM — galaxies: starburst — ISM: 



1. INTRODUCTION 

Our understanding of the cosmic star formation history 
(SFH) of gal axies has progr e ssed significantly over the past 
decade (e.g., Hopkins 2004; Hopkins & Beacom 2006). In 
the same time the space density of neutral hydrogen gas has 
been mea sured over the majority of cosmic history (see Fig- 
ure 8 of|Lah et al.ll2007l) . The evolution of the atomic hydro- 
gen (HI) in the universe will be comprehensively determined 
within the next few years by extremely sensitive surveys with 
the next generation of radio telescope instrumentation (e.g., 
van der Hulst et al.ll2004T: iRawlings et al.ll2004t Uohnston etafl 
2008), and it is timely to consider mechanisms associated with 
this evolution. 

The space density of HI in galaxie s appears to evo lve sur- 
prisingly little from z « 5 to z rs 0.2 dLah et al.ll2007l) . a span 
of roughly lOGyr, the latter half of which sees a decline in the 
space density of star formati on rate (SFR) in galaxies b y almost 
an order of magnitude (e.g., Hop kins & Beacorrj|2006l) . Given 
the SFR density it is easy to show that the HI plus molecu- 
lar gas at high redshift would be exhausted on timesc ales of 
a few Gyr if it were not continually replenished. lErbl ([2008) 
presents a model incorporating gas infall, outflows and con- 
sumption by star formation, to explain both replenishment and 
the mass-metallicity relation in high-redshift (z rs 2) galaxies. 

Hydrodynamic simulations advocating hot and cold modes 
of accretion indicate that the infall rat e closely tracks the sta r 
formation rate (e.g. iKeres et ail 120051: iBirnboim et al.l 12007). 
with star formation moderated by the rate of infall. The sim- 
ulations, however, neglect gas outflows from galaxies, which 
are a significant component of gas depletion. The quantitative 
infall rates predicted are thus insufficient to maintain a constant 
HI density in galaxies. Observed rates of gas infall in local 
galaxies, also, are only about 10% of the star formation rate 
(San cisi et al.ll2008l) . The difficulties in explaining replenish- 
ment through infall leave the physical mechanism of this re- 
plenishment as a critical open question in galaxy evolution. 

In this Letter we suggest a mechanism directly associated 

1 Thanks to Sandhya Rao (Rao et al. 2006) for this terminology. 



with the SFR in galaxies that can provide the necessary replen- 
ishment of neutral gas to maintain an essentially unevolving, or 
slowly evolving, HI mass density. We infer the density of gas 
required to reproduce the observed SFH in §|2] In §|3]we present 
a number of models for the replenishment of this gas, and show 
that a replenishment proportional to the SFR density can repro- 
duce the necessary gas mass density. A replenishment mech- 
anism associated with galactic supershells is detailed in §[4] 
and the results are summarised in §0 Throughout this analysis 
we adopt the "737'Q cosmol ogy with Hp = 70 kms~ 1 Mpc -1 , 
tt M = 0.3, n A = 0.7 (e.g. JSpergel et alj|2003l) . 

2. ESTIMATING THE MASS DENSITY OF STAR FORMING GAS 

While our motivation is to understand the observed lack of 
significant evolution in the HI mass density, we approach this 
by considering the total mass density of gas available to form 
stars, which includes molecular as well as atomic gas. We ne- 
glect the intricacies in the conversion of HI to molecular gas 
associated with the star formation process, as this occurs on 
timescales very short compared to those involved in this analy- 
sis. What is important is the total reservoir of gas available for 
star formation, "star-forming gas," psfg, at a given redshift, 
comprised of both atomic and molecular gas. Star formation 
is an inefficient process and psfg will likely be somewhat less 
than the total of the atomic and molecular gas mass densities. 

The cosmic evolution of the HI mass density has been chal- 
lenging to measure. There is an implicit assumption that the 
damped Lya absorbers used to measure the HI mass density at 
z > 0.3 are representative of the galaxies used to trace the SFH. 
The best current measu rements suggest that this i s not an un- 
reasonable assumption CZwaan & Prochaskall2006l) . while also 
highlighting the difficulty in obtaining observational constraints 
on psfg- In the absence of direct observational measurements, 
we estimate psfg indirectly from the observed SFH, using the 
local relat i onship between gas and SFR surface densities from 
Kennicutt dl998l). W e calculate psfg using Equation (10) of 
Hop kins et al.1 (120051) . and the results are shown in Figure Q] 
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The open triangles in Figure Q] ar e derived from the piecewise 
linear fit to the SFH of Hopkins & Beacoml <HQ06), their Ta- 
ble 2, for their "SalA" initial ma ss function (IMFfl T he open 
squares correspond to the SFH of Wil kins et alj d2008l) . This is 
based on an evolving IMF, more top-heavy at higher redshift, 
which reduces the derived SFR densities and gives lower in- 
ferred psfg at high-z than a universal IMF. 

These estimates make the implicit a ssumption that the 
Kennicutt-Schmidt law for star formation ( iKennicuttll 19981) is 
valid at all re dshifts, an assumption questioned by a num- 
ber of authors dZwaan & Prochaskal2006tlWolfe & Cherj2 006; 
IWild et al]|2007l) . The relation has, though, been successfully 
used to expl ain t h e mas s-metallicity relationship seen in z « 2 
galaxies by lErbl d2008l) and in hydrodynamic simulations by 
Rob ertson & Kravtsovl d2008l) seems to arise as a natural con- 
sequence of gas physics. In the absence of detailed observations 
of molecular gas in large numbers of high redshift galaxies, the 
assumption that the local relationship holds at higher redshift is 
the simplest with which to explore further. 

The gas consumption timescale comes from dividing psfg 
by the stellar mass density remaining in stars and stellar rem- 
nants, inferred from the SFR density. This gives timescales of 
w 1-5 Gyr with the lower value at higher redshift where the 
SFR density is largest. This short timescale, derived from the 
global properties of all galaxies, is strikingly similar to the gas 
consumption timescales fo und for individual nearby galaxies 
dWong & Blitdl200Tl [2002b . 

3. THE EVOLUTION OF THE STAR FORMING GAS 

To explore the interaction between gas consumption and 
replenishment we show the effects of several simple mod- 
els in Figure [2] All models begin at a lookback time, 
£l, of 12.55 Gyr (z = 6) assuming an initial value for 
log(psFG/(M Mpc~ 3 )) = 8.1, with Psfg(^l) calculated sim- 
ply as the integral over time. At each time-step gas is con- 
sumed by star formation (— p*(t)) and a similar amount in gas 
outflows^. Gas is returned to the interstellar medium (ISM) 
through stellar evolutionary processes (stellar wi nds, supernova 
ejecta ) with a recycling fraction of R (e.g ., iKennicutt et alj 
Il994t iMadau et al.1 119981: ICole et alj l200lh . adding a factor 
+Rp*(t). For the SalA IMF, R = 0. 4 while other IMFs 
will h ave different recycled fractions (e.g., Hop kins & Bea com 
2006). Finally a replenishment factor K(t) is added. This can 
be expressed as: 



and K(t) = 0.21 M© Mpc -3 yr" 1 (dash-triple-dotted line), are 
also shown. The "bounce" seen in these models arises from an 
early excess in consumption followed by progressively decreas- 
ing consumption as declines for z < 1, predicting excess 
Psfg at lookback times t L < 4 - 6 Gyr (z = 0.4 - 0.7). 

Replenishment factors proportional to the SFR density pro- 
vide an obvious way to ensure psfg remains constant with 
time. A replenishment factor of K (t) — 1.6/9, in Equation (fl~|i 
gives ApsFG = at all redshifts (dotted line in Figure [2). 
Different constants of proportionality allow for slowly varying 
changes in psfg- With K(t) = 1.5 2/7* (a fact or 0.95 of that 
required for complete replenishment; Erb 2008) the replenish- 
ment does not fully balance consumption, and gives a slow de- 
cline in the global gas density (heavy solid line in Figure[2]i. 

In the following section we consider supershells in the ISM 
as the driver of a physical mechanism for replenishment. Super- 
shells are both directly associated with star formation and are 
highly efficient at converting hot phase gas into star-forming 
gas. 

4. DISCUSSION 

4.1. A physical mechanism for replenishment 

The ISM in galaxies contains expanding supershells or super- 
bubbles associated with previous generations of star formation. 
We propose that the neutral and molecular gas replenishment 
in the walls of supershells is sufficient, and of the appropriate 
form, to provide a natural mechanism explaining the relatively 
flat evolution in the HI mass density. 

Supershells are large scale expanding shells of gas driven 
by s upernovae (SNe) and stellar winds from OB star clus- 
ters dOevlll9"96l:IOev & Smedlevlll998tlMcClure-Grifnths et al.1 
2001). Supershells have long been suggested to have a trig- 
gering effect on subsequent generations o f star formation 
(Mc Crav & Kafatos) 119871: lElmegreenl 119981: lHartmann etafl 



2001; Bergin et alj (2004t lOev et all 120051) . and are effective 



Psfg(*l) = PsFc(i = 12.55) + 

rt=th 

(-1.6p*(t) + K(t))dt. (1) 

/t=12.55 

These simplifications hide a wealth of complex ISM and 
IGM interactions, including the fact that material returned to the 
ISM through stellar evolution, as well as infalling gas, may con- 
tribute to ionized, or otherwise non-star-forming components as 
well as directly to psfg- These details are subsumed into the 
effective replenishment factor, K(t). 

A model with no gas replenishment, K (t) — 0, is shown 
in Figure [2] (the dashed line), emphasising the rapid consump- 
tion timescale. Two models assuming constant rates of re- 
plenishment, K(t) = O.16M Mpc~ 3 yr _1 (dash-dotted line) 

2 The SalA IMF is a modified Salpeter IMF with a turnover below 0.5 M Q (detailed in lHopkins & Beacoml200r3 and lBaldrv & Glazebrookll2003l) . 

3 Wefollow lErbl<200a) in assuming the gas outflow rate from winds exhausts an amount equal to the SFR, i.e., — p„ (t) . This is a coarse approximation, but based on 
galactic outflow rates in the nearby universe I Veilleux et al. 2005) it may be reasonable when averaged over large populations of star forming galaxies. 

4 Here and throughout we assume the inclusion of all core-collapse supernovae: types II, lb, and Ic 



at replenishing star-forming gas through several mechanisms. 
First, supershells are efficient at cooling and recombining ion- 
ized gas through radiative cooling in shell walls dKoo et ail 
[1991 . This may be critical in converting gas shock-heated 
by previous generations of SNe within a galaxy, or new hot- 
mode infall gas, to a potentially star-forming state, as the 10 6 K 
gas may otherwise never cool to support subsequent star for- 
mation. Second, molecular gas can form from neutral gas in 
situ in shell walls, where compression and the development of 
instabilities leads to sufficiently high neutral gas densities to al- 
low for cooling and self-shielding on timesc ales of 10 — 20 Myr 
dBergin et alj|2004t iHennebelle et al.l l2008). Finally, they can 
compress pre-existing molecular mat erial to trigger m olecu- 
lar cloud collapse and star formation ( El megreenll 1 9981) . The 
timescales for these processes are shorter than, or comparable 
to the supershell lifetime (~ 20 Myr), which is in turn short 
compared to the global gas consumption timescale of several 
Gyr. 

To establish whether the replenishment achievable in su- 
pershells is sufficient to make this mechanism feasible, we 
first convert the replenishment rates given in §[3] into a replen- 
ished mass per SN event. A replenishment rate proportional 
to the SFR density is also proportional to the rate of super- 
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nova type II (SNlfl). Converting a proportionality to SFR den- 
sity into one depending on the SNII rate, <6sntt, d epends on 
the assumed IMF. From Hopkins & Beacom (2006) psnii = 
(O.OO915/M )p» for the SalA IMF. The replenishment rate 
K(t) = 1.6/6, becomes K(t) = 174.9/6 SNn M . 

The other extreme ch oice of IMF consistent wi th the nor- 
malization of the SFH dHopkins & Beacoml 1200 6) is that of 
Bald rv & Glazebrookl d2003l hereafter the BG IMF). For the 
BG IMF p snii = (0.0132/Mq) p, The recycled fraction is 
R = 0.56 (Hopkins & Beacom 2006), changing the consump- 
tion term in Equation (Q~|i to —1.44/6,. The corresponding re- 
plenishment rate is K (t) — 109.1/6snii Mq. These extremes 
imply that, depending on the IMF, sufficient gas replenishment 
to maintain a constant HI mass density with redshift would be 
achieved if each SN event caused the recombination and cool- 
ing of s» 110 — 180 Mq of gas. These IMFs are the extrema 
given the SFH normalization limits, and most reasonable IMFs 
should result in masses within this range. 

Detailed measurements to confirm molecular gas formation 
within supershells are observationally challenging. We use 
the limited data currently available to assess the replenishment 
rates associated with supershells, and to establish whether at 
le ast one well-studied supers hell a chieves the req u ired ra te. 

iMcClure-Griffithsl d2005l) and iDawson et all d2008l) have 
shown explicit cases of molecular clumps along the edges of su- 
pershells, suggestive of some degree of in situ formation, with 
a significant amount of molecular material asso ciated with the 
supers hell walls. The supershell investigated bv lDawson et al.l 
(2008) is associated with about 2 x 10 5 Mq of molecular gas, of 
which those authors estimate that 80% likely comes from a pre- 
existing giant molecular cloud. Of the remaining « 4 x 10 4 Mq 
of molecular gas it is difficult to determine how much is pre- 
existing and how much has been cooled and recombined by 
the expansion of the shell. We can use w 4 x 10 4 Mq as 
an upper limit to the replenishment rate. About 30 stars with 
stellar mass M* > 7 Mq are required to form this supershell, 
including stars that may not yet have gone supernova. This 
gives < 1300 — 2000 Mq of molecular mass replenished per 
SN event, a limit comfortably encompassing the required rate. 
This upper limit could change significantly depending on the 
fraction of pre-existing molecular material and also on the frac- 
tion of stars that have not yet gone supernova. 

Not all SNe l ie within supershells, although 

Hig don & Lingenfelterjd2005l) estimate that a minimum of 65% 
of SNII should occur in superbubbles, increasing to ~ 80 — 90% 
when the spatial and temporal correlations of stellar clusters are 
considered. If 80% of SNII are associated with supershells, for 
example, this would increase the required replenishment rate 
per SN to w 140 - 230 Mq. But even if as few as 10% of 
all SNII contribute in thi s way to the r e plenis hment, the rate 
implied by the results of IDawson et al.l d2008l) would still be 
sufficient. This confirms that the necessary replenishment rates 
are likely to be achievable within supershells. 

The observed decline by a factor of two in the HI mass 
density may be a natural consequence of a replenishment rate 
about 95% of that required to match consumption, as shown 
by the heavy solid line in Figure [2] If the actual replenishment 
rate from supershells lies somewhere between the required rate 
and our derived upper limit, though, there may in fact be too 
much newly replenished gas to allow any decline in the neu- 
tral gas mass density. A possible resolution in this scenario 
would be increasing the proportionality between the gas out- 



flow rates and the SFR as redshift decreases. This is not unrea- 
sonable, as the SFH is becoming progressively mo re dominated 
by lower-mass galaxies with decreasing redsh ift dJuneau et al.l 
120051; iPanter et alj|2007l: iMobasher et alj|2008l) . Galaxies with 
stellar masses M, < 10 10 M Q dominate the SFH at z < 1 
(Mobashe retail 120081) . Such low-mass galaxies lose more 
mass in gas outflows in proportion to their SFR than high- 
mass galaxies, simply due to the former's shallower poten- 
tial wells (e.g.,|5ek el & Silkl ll986t iMac Low & Ferraralll999 l : 
iFerrara & Tolstovl2 000). This effect may contribute to the slow 
decline in the HI mass density. 

4.2. Limitations of the proposed mechanism 

We have treated a number of complex physical processes in 
very general terms. While being cautious of oversimplification, 
we have attempted to capture the essential interactions between 
star formation, recycling from stellar evolutionary processes, 
ISM processes of heating and ionization, recombination, cool- 
ing and molecule formation, together with infall from the IGM, 
and outflow of ISM material. Most of this complexity is con- 
cealed within the replenishment factor K(t). 

One issue is that stellar winds and SNe contribute to all com- 
ponents o f the ISM rather than solely to psfg- In a "gal actic 
fountain" dShapiro & Field! 1976tlHouck & Bregmanl 19901) . in- 
falling gas will contribute to, and outflowing gas will strip from, 
all components. If recycled gas includes a component that never 
subsequently forms stars (such as some recycled gas in the ion- 
ized phase being ejected from the galaxy before contributing 
to star formation), the factor +i?p* it) in Equation (Q]) will be 
reduced, and K(t) will need to be increased to compensate. 

Our quantitative results strongly depend on the assumed gas 
outflow rate. Variations by a factor of two or so in either di- 
rection will still result in a constant or slowly varying HI mass 
density, as long as a proportionalit y with the SFR of th e host 
galaxies remains (as suggested by IVeilleux et alj|2005l) . The 
chosen outflow rate is an effective average over all star forming 
galax i es, and is consistent with observed tre nds (e.g., iMartinl 
119991: iPettini et al.|[2000i: IVeilleux et al1l2005l) . While individ- 
ual galaxies show a large observed scatter between outflow rates 
and SFRs, for the ensemble properties of the total population 
this assumption should be robust. 

The proposed replenishment through the supershell mecha- 
nism is not inconsistent with some simultaneous replenishment 
through infall. Metallicity considerations, which we do no t ad- 
dress here, do require infall of some low metallicity gas jErbj 
2008l j), and gas infall in loca l galaxies is well est ablished (e.g., 
Bland-H awthorn et alj|2007l: ISancisi et al.ll2008l) . although the 
observed infall rate is insufficient to match consumption. 

5. SUMMARY 

We have identified a possible resolution to the puzzle of why 
the HI mass density of the universe evolves so little for so much 
of cosmic history. We propose that replenishment is driven by 
supershells associated with star forming complexes in galax- 
ies. Pre-existing ionized gas efficiently cools and recombines 
in supershell walls. Molecular gas forms in situ in shell walls, 
and molecular material is compressed to trigger cloud collapse 
and star formation. This mechanism provides a natural expla- 
nation for replenishment that has the desired proportionality to 
the SN rate . The level of replen ishment observed in a Galactic 
supershell dDawson et al . 2008) appears more than sufficient to 
provide the required replenishment rate of « 110 — 180 Mq per 
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SN event. 

The factor of two decline in the HI density between z«0.2 
and z = could be explained through either (1) a replenish- 
ment rate that is marginally lower than that required to ex- 
actly balance gas consumption; or (2) the inability of low- 
mass galaxies, which dominate the star formation history in this 
epoch, to retain their newly formed HI; or, perhaps more likely, 
a combination of both. 



We thank the referee for their comments and helpful input 
on hydrodynamic simulations. We thank Joss Bland-Hawthorn 
and Leo Blitz for valuable discussions that helped to improve 
this analysis. AMH acknowledges support provided by the 
Australian Research Council (ARC) in the form of a QEII Fel- 
lowship (DP0557850). BMG acknowledges support from the 
ARC through a Federation Fellowship (FF0561298). 



REFERENCES 



Baldry, I. K., Glazebrook, K. 2003, ApJ, 593, 258 

Bergin, E. A., Hartmann, L. W., Raymond, J. C, Ballesteros-Paredes, J. 2004, 
ApJ, 612, 921 

Birnboim, Y., Dekel, A., Neistein, E. 2007, MNRAS380, 339 
Bland-Hawthorn, J., Sutherland, R., Agertz, O., Moore, B. 2007, ApJ, 670, 
L109 

Cole, S., et al. 2001, MNRAS, 326, 255 

Dawson, J. R., Mizuno, N., Onishi, T., McClure-Griffiths, N. M., Fukui, Y. 

2008, MNRAS, (in press; arXiv:()802.4463 I 
Dekel, A., Silk, J. 1986, ApJ, 303, 39 

Elmegreen B. G., 1998, in "Origins," eds., Woodward C. E., Shull J. M., 

Thronson Jr. H. A., ASP Conf. Ser. 148, 150 
Erb, D. K. 2008, ApJ, 674, 151 
Ferrara, A., Tolstoy, E. 2000, MNRAS, 313, 291 

Hartmann, L., Ballesteros-Paredes, J., Bergin, E. A. 2001, ApJ, 562, 852 
Hennebelle, P., Banerjee, R., Vazquez-Semadeni, E., Klessen, R., Audit, E. 

2008, A&A, (submitted; arXiv:0805. 13661 
Higdon, J. C, Lingenfelter, R. E. 2005, ApJ, 2005, 628, 738 
Hopkins, A. M. 2004, ApJ, 615, 209 
Hopkins, A. M., Beacom, J. F. 2006, ApJ, 651, 142 
Hopkins, A. M., Rao, S. M., Turnshek, D. A. 2005, ApJ, 630, 108 
Houck, J. C, Bregman, J. N. 1990, ApJ, 352, 506 
Johnston, S., et al., 2008, PASA, 24, 174 
Juneau, S., et al. 2005, ApJ, 619, L135 

Keres, D., Katz, N., Weinberg, D. H., Dave, R. 2005, MNRAS, 363, 2 
Kennicutt, R. C, Jr. 1998, ApJ, 498, 541 

Kennicutt, R. C, Jr., Tamblyn, P., Congdon, C. E. 1994, ApJ, 435, 22 
Koo, B.-C, McKee, C. F. 1992, ApJ, 388, 93 
Lah, P., et al. 2007, MNRAS, 376, 1357 

Madau, P., della Valle, M., Panagia, N. 1998, MNRAS, 297, L17 

Mac Low, M.-M., Ferrara, A. 1999, ApJ, 513, 142 

Martin, C. L. 1999, ApJ, 513, 156 

McCray, R„ Kafatos, M. 1987, ApJ, 317, 190 



McClure-Griffiths, N. M. 2005, in "New Horizons in Astronomy," eds S. J. 
Kannappan, S. Redfield, J. E. Kessler-Silacci, M. Landriau, N. Drory, ASP 
conf. ser. 352, 95 

McClure-Griffiths, N. M., Dickey, J. M., Gaensler, B. M., Green, A. J. 2001, 

ApJ, 562, 424 
Mobasher, B. t et al., 2008, ApJ, (submitted) 

Oey, M. S., Watson, A. M., Kern, K., Walth, G. L. 2005, AJ, 129, 393 

Oey, M. S. 1996, ApJ, 465, 231 

Oey, M. S., Smedley, S. A. 1998, AJ, 116, 1263 

Panter, B., Jimenez, R., Heavens, A. E, Chariot, S. 2007, MNRAS, 378, 1550 
Pettini, M., Steidel, C. C, Adelberger, K. L., Dickinson, M., Giavalisco, M. 

2000, ApJ, 528, 96 
Prochaska, J. X., Herbert-Fort, S., Wolfe, A. M. 2005, ApJ, 635, 123 
Rao, S. M., Turnshek, D. A., Nestor, D. B. 2006, ApJ, 636, 610 
Rawlings, S., Abdalla, F. B., Bridle, S. L., Blake, C. A., Baugh, C. M., 

Greenhill, L. J., van der Hulst, J. M. 2004, New Ast. Rev., 48, 1013 
Robertson, B. E., Kravtsov, A. V. 2008, ApJ, (submitted; arXiv:0710.2102 1 
Sancisi, R., Fraternali, F, Oosterloo, T, van der Hulst, T. 2008, A&A Rev., (in 

press; arXiv:0803.0109 1 
Shapiro, P. R., Field, G. B. 1976, ApJ, 205, 762 
Spergel, D. N. et al., 2003, ApJS, 148, 175 

van der Hulst, J. M., Sadler, E. M., Jackson, C. A., Hunt, L. K., Verheijen, M., 

van Gorkom, J. H. 2004, New Ast. Rev., 48, 1221 
Veilleux, S., Cecil, G, Bland-Hawthorn, J. 2005, ARA&A, 43, 769 
Wild, V., Hewett, P. C, Pettini, M. 2007, MNRAS, 374, 292 
Wilkins, S. M., Trentham, N., Hopkins, A. M. 2008, MNRAS, 385, 687 
Wolfe, A. M., Chen, H.-W. 2006, ApJ, 652, 981 
Wong, T., Blitz, L. 2001, Ap&SS, 276, 659 
Wong, T., Blitz, L. 2002, ApJ, 569, 157 
Zwaan, M. A., Prochaska, J. X. 2006, ApJ, 643, 675 

Zwaan, M. A., Meyer, M. J., Staveley-Smith, L., Webster, R. L. 2005, MNRAS, 
359, L30 



Evolution of gas and star formation 



0.2 



0.5 



2 3 



x 

02 



O 



o o 



- 

m 
Q. 

^ o 

X 








lookback time (Gyr) 



10 



FIG. 1 . — The cosmic history of neutral and "star-formin g" gas mass density. Solid circles : The HI density as shown in Figure 8 of Lah et al. (2007). From high 
to low redshift the data come from Proch aska et alj {2005) at t L > 10 Gyr (z > 1.9); from lRao et"aT] <2006D between 5 < t L < 9 G yr, (0.61 < z < 1 2 2); the 
stacking measurement from Lah et al. (2007j) at 4l = 2.8 Gyr, (z = 0.24); and the HIPASS m easurement of the l ocal HI density from Zwaan et al. ( 2005). Open 
red triangles: PsfG> the neutral plus molecular gas density inferred from the SFH from Hopkins & Beacom 12006), assuming the Kennicutt-Schmidt relation for 
star formation; Open blue squares: psfg as inferred from the SFH of Wilkins et al. 1 2008j), which assumes an evolving initial mass function. The solid magenta 
and cyan lines are nominal parameterisations of the possible evolution of the total gas reservoir for each of the SFH cases. 
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FIG. 2. — Models for the evolution in the gas reservoir. The HI mass density (solid circles), psFG (open blue squares), and solid cyan line are as in Figure[T] 
Dashed line: Predicted gas reservoir evolution assuming no replenishment; Dash-dot and dash-triple-dot lines: Two different rates of constant replenishment; Dotted 
line: Replenishment rate proportional to the SN rate; Heavy solid line: Replenishment rate proportional to the SN rate, but at a factor of 0.95 of that required to 
balance the consumption. See text for further details. 



